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in  a  local  verticE.l  coordinate  system,  two  sets  of  elevation  angles  are 
found  for  the  tumbling  rocket,  and  the  set  that  satisfies  an  initial  condition 
exactly  is  the  correct  one.  The  deviation  of  the  rocket  tumbling  plane  from 
the  vertical  plane  can  be  determined  at  the  times  nf  closest  encounter  of  the 
magnetic  field  vector  and  rocket  vector.  The  deviation  at  all  times  can  be 
approximately  determined  b>  si  i.tsoidal  curve  fitting  with  tumbling  frequency. 
The  angles  hetween  the  rocket  vector  and  moonlight  vector  during  the  instants 
of  moonviews  at  the  onboard  camera  are  also  calculated. 
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Attitude  Determination  of  a  Spinning  and  Tumbling  Rocket 
Using  Data  From  Two  Orthogonal  Magnetometers 


1.  INTRODUCTION 


This  work  was  motivated  by  data  obtained  on  Air  Force  scientific  rocket  No. 

A3 1.603.  Electron  and  ion  beam  experiments1  were  carried  out  on  the  rocket.  In 
order  to  analyse  and  interpret  the  scientific  data,  it  is  necessary  to  know  the  atti¬ 
tude  of  the  rocket  during  the  instants  of  beam  emissions.  Of  particular  importance 
is  the  angle  between  the  rocket  and  geomagnetic  field  lines,  because  the  presence 
of  the  geomagnetic  field  gives  rise  to  anisotropic  electron  motion.  The  electrons 

move  easily  along,  but  with  difficulty  across  the  geomagnetic  field.  Furthermore. 
2  3 

there  are  indications  ’  that  the  magnitude  of  return  current  is  less  for  a  beam 
emitted  downwards  along  the  magnetic  field. 

Other  interests  also  motivated  this  work.  The  standard  technique  of  attitude 
determination  with  magnetometer  data  requires  data  from  three  orthogonal 

(Received  for  publication  31  March  1981) 

1.  Cohen.  H.  A.,  Sherman.  C. ,  and  Mullen.  G.  (1979)  Spacecraft  charging  due  to 

positive  ion  emission:  an  experimental  study,  Geophys.  Res.  Lett. 
jJJNo.  8):515.  -  - 

2.  Winckler.  J.R.  (1980)  The  application  of  artificial  electron  beams  to  magneto- 

spheric  research.  Review  of  Geophysics  and  Space  Physics,  l£:659-682. 

3.  Jacobsen,  T.A. ,  and  Maynard,  N. C.  (1980)  Polar  5— an  electron  accelerator 

experiment  within  an  aurora  3.  Evidence  for  significant  spacecraft  charging 
by  an  electron  accelerator  at  ionospheric  altitudes.  Planet.  Sp.  Sci. 

££29 1-307.  - - 
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magnetometers,  the  magnetometer  telemetry  signal  conversion  factors  being  known 
from  pre -launch  calibration.  However,  when  only  data  Trom  two  orthogonal  mag¬ 
netometers  are  available  and  the  conversion  factors  are  not  given,  it  becomes  a 
challenging  problem  whether  the  attitude  can  still  be  estimated  in  a  local  vertical 
coordinate  system.  For  an  object  tumbling  slowly,  compared  with  its  spinning,  in 
a  vertical  plane  with  small  sinusoidal  variations,  a  method  is  offered  for  solving 
this  problem. 

The  position  £  of  the  spacecraft  and  the  Earth's  magnetic  field  vector  B  at  that 
position  are  given.  The  magnetometer  data  (hereafter  abbreviated  ad  m.d. )  are 
known  to  be  approximately  linearly  related  to  the  true  magnetic  data  B  by: 

(M  -  D)S  =  B 

where  M,  JD»  and  S  are  m.  d. ,  bias  offset,  and  scaling  factor  respectively.  For  a 
triaxial  magnetometer  set,  the  three  components  of  M  can  be  obtained.  If  only  two 
magnetometers  are  available,  additional  techniques  have  to  be  devised.  The  main 
objective  is  to  determine  the  directions  of  the  magnetometers  with  respect  to  the 
local  vertical  coordinate  system,  those  with  respect  to  the  spacecraft  coordinate 
system  being  given. 


2.  BIAS  OFFSET 

If  the  m.d.  varies  in  a  simple  Sinusoidal  manner  and  then  for  a  sufficiently 
long  period  in  a  steady  state,  the  bias  offset  can  be  determined  as  the  mean  of  the 
m.d.  in  that  period,  because 

<  A  sin  w(t  - 1  )  +  D>  =  <  A  sin  (t  -  X  )>+D*D. 
o  o 

If  more  than  one  significant  frequency  is  present  in  the  steady  state  time  series, 
the  bias  offset  can  be  determined  in  the  same  manner,  because 

Lim  <  Zy  A.  sin  u*U  -t)  +  D>xS<A.  sin  ui.(t  -  t  )  >  +  D  «  D  , 

t  "**  i  1  ^  °  i  1  1  ° 

This  method  applies  to  the  time  series  of  each  magnetometer. 

If  the  time  span  of  available  data  is  short,  such  as  three  periods  (3  X  a»”S  or 
less,  for  example,  then  the  offset  D  has  to  be  determined  by  least  square  curve 
fitting,  if  good  accuracy  is  desired.  For  rigid  bodies,  sine  functions  with  known 
frequencies  are  usually  suitable  for  fitting  purposes.  The  frequencies  can  be  com¬ 
puted  by  means  of  fast  Fourier  transform  or  maximum  entropy  methods. 
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S.  SCAI4NG  FACTOR 


f  1 


After  removing  the  bias  offset,  D,  the  magnetometer  data  {Ml  become  direct¬ 
ly  proportional  to  the  magnetic  field  components  {  b}  .  If  data  from  only  two  mag¬ 
netometers  were  available,  it  would  not  be  possible  to  derive  exact  knowledge  of 
the  three  dimensional  vector  M  ft)  at  any  time  t.  However,  in  an  ideal  situation— that 
is,  when  the  spin  vector  of  the  vehicle  remains  constant  in  direction  with  respect 
to  an  inertial  coordinate  system  in  which  the  magnetic  field  lines  are  at  rest— the 
scaling  factor  can  be  determined  as  follows. 

Let  Mjft)  and  Mj(t)  be  available  during  a  spin  period.  Mg(t)  being  unavailable. 
Let  Mjit)  be  the  measurement  in  the  spin  direction.  The  Earth's  magnetic  field 
B(t)  is  assumed  constant  in  the  space  and  period  considered.  With  the  vehicle 
spinning  about  the  one -direction,  the  magnetometer  readings  }M(t)f  at  time  t  are 
given  by 


lM(t>]  =10  0 

0  cos  p(t)  sin  p(t) 

0  -sin  p(t)  cos  9(t) 


lM(t0>] 


where 


|  B|  =  S  |  NlJ(to>  +  Mj(tQ)  +  M*ft0>J.J 


S  being  the  scaling  factor. 

In  a  spin  period,  the  extremum  value  of  Mjit)  is  given  by 
8Ma(t)/»d(t)  *  0 


which,  using  Eq.  (1),  gives 

9<t)  =  tan"1  £  M3<to)/M2<Vj 


Max[Ma<t>]  =  |Ma(t0>  +  M^iJ1 


so  that 


t?(t)  +  Max2 


M}‘/2  ■  “2«„>  *  “j«o'}U! 


Therefore* 


|b| 


{“! 


(t)  ♦  Mm*  (Ma<t» 


(4) 


Thus,  despite  the  lacking  of  Mj't)  dete  et  any  time  t.  the  seeling  lector  S  can  be 
determined,  using  fiq.  (4).  if  the  value  1  B|  is  known  by  using  standard  model  cal¬ 
culations  or  other  appropriate  data  if  available. 


4.  ASPECT  ANGLE 


The  aspect  angle  a  ^t)  subtended  by  the  ith  magnetometer  vector  with  the 
Earth’s  magnetic  field  line  B  at  time  t  is  given  by 


a  jit) 


cos 


-1 


Thus. 


a.(t)  «  cos 


-1 


SM^t)*  B(t) 
I  Bit)  |* 


M^t) 


(5) 


VwJ(t)  ♦  Max2  |^M2(tl 


(6) 


Of  particular  Interest  for  Rocket  A31.603  is  magnetometer  No.  1,  which  points 
forward  along  the  rocket  axial  vector,  the  spin  axis  being  about  coincidental  with 
the  geometric  axis.  Thus,  the  aspect  angle  O  ^{X)  gives  the  angle  subtended  by  the 
rocket  axis  with  the  B  line  at  the  rocket's  location. 


5.  TUMBLING  AND  SPINNING  ROCKET 


In  general  situations,  the  spin  vector  of  a  vehicle  does  not  remain  constant  in 
direction  with  respect  to  an  inertial  coordinate  system.  For  a  pure  tumbling  and 
spinning  vehicle,  the  magnetometer  data  { M(t) }  measured  on  the  vehicle  at  time 
t  are  given  by 


[M(t>] 


’  1 

0 

0 

cos  8  (t) 

0 

sin 

8(t) 

0 

cos  p{t) 

sin  p(t) 

0 

1 

0 

0 

-sin  p(tl 

cos  p(t) 

-sin  9  (t) 

0 

cos 

9(t) 

[M(t=0)J 


(7) 


4.  SUA/AFGL  (1979)  Data  Report  Rocket  No.  A31.  603. 
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where  the  tumbling  axis  la  the  y-axls  In  the  Inertial  system  end  the  spinning  eels 
ts  the  x«axia  In  the  spacecraft  system.  At  time  t  •  0,  9  (t)  end  p(t)  ere  defined 
as  0.  In  a  spinning  period,  the  maximum  value  of  Mg(t)  is  measured  when  p(t) 
satisfies  the  following  equation,  at  t  *  tm« 

£ -sin  9  sin  f  Mj(0)  4  cos  p  Mj(0)  4  cos  9  sin  $  M3(0)J  •  0 


which  gives,  using  Eq.  (7). 

Max  [Ma(t>]  «  |  m|(0>  4  cos2  «<tm)  mJ(0>  4  sin*  0  <tm> 

3  >1/2 

Mf(0>  -  2  cos  «(tm)  sin  #(tm>  M^O)  Mj(O)  S  .  (8) 

If  the  spin  frequency  is  much  higher  than  the  tumble  frequency,  then  9  (tm>  would 
not  be  too  different  from  9  (t*  0).  which  is  r.ero.  Thus, 

Una  Max[Ma<*>]  -  j  M*(0)  4  mJ<0)|1/2  ,  (9) 

which  reduces  to  Eq.  (3).  the  result  for  no  tumble  case. 


6.  ERROR  ANALYSIS  ON  ASPECT  ANGLE  ALGORITHM 

For  low  tumble  frequency /spin  frequency  ratio,  the  aspect  angle  a  .(t)  sub¬ 
tended  by  the  ith  magnetometer  vector  with  the  B  line  is  given  by,  using  Eq.  (6), 


where  MaxlMgftll  is  given  in  Eq.  (8).  This  algorithm  has  its  limitation,  and  will 
now  be  discussed.  The  term  Max[Ma(t)]  is  used  in  a  .  formula  because  of  necessity 
due  to  the  unavailability  of  Mg(t)  data.  If  there  is  no  spin  vector  movement,  it  has 
been  shown  (in  Eq.  (3)]  that  this  term  really  equals  the  vector  sum  of  the  B  compo¬ 
nents  perpendicular  to  the  spin  axis.  However,  if  the  spin  axis  moves,  then,  as 
shown  in  Eq,  (8),  the  term  Max(Mg(t)]  becomes  a  not  too  simple  combination  of 
M2(tc).  M3(tQ).  and  even  M^).  The  error  in  a^(t)  is  given  by 
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6,a  (t)  =  cos 


-1 


M^t) 


V 


mJ(0)  +  M2(0)  +  M2(0) 


“COS 


M^t) 


V 


M?(t  )+  M?(0)+  cos2  Bit) 

l  zn  &  m 


M?(0)  +  sin2  Bit)  M?(0)  -  2  cos  Bit)  sin  Bit)  M,(0)  M,(0) 
3  mi  m  m  l  s 


(10) 


which  is  time  (t  and  t  )  dependent.  For  given  values  of  M.(t),  M,(0),  and  6  (t  ),  the 
m  11  m 

error  Aa  ^(t)  can  be  computed,  using  Eq.  (10).  For  a  time  span  longer  than  a  spin 

period,  it  is  reasonable  to  expect  that  &a  (t)  is  distributed  with  its  mean  near  zero. 

For  Rocket  A31.603,  the  value  of  <6( t_)>  *  18®.  This  can  be  observed  from  the 

m 

variations  of  a  ^(t)  and  a  2(t):  in  a  period  of  a  2(t),  which  is  the  aspect  angle  of  the 
magnetometer  in  spin  plane,  there  are  approximately  18®  variations  of  a  ^ (t),  which 
is  the  aspect  angle  of  the  magnetometer  looking  forward  along  the  spin  axis.  If  the 
values  of  all  the  terms  in  Eq.  (10)  are  known,  the  error  Aa.(t)  can  be  determined 
at  any  time  t. 

The  realistic  situation  is  more  complicated,  because  the  vehicle  not  only 
tumbles  but  it  precesses  in  a  large  conic  angle.  In  general,  three  Euler  angles 
are  needed  for  arbitrary  rotations.  Accordingly,  Eq.  (7)  has  to  be  modified  to 
accommodate  three  rotation  matrices  on  the  right-hand  side,  and  Eq.  (10)  becomes 
much  more  complicated.  If  such  a  detailed  analysis  is  necessary,  it  can  be  done. 
However,  if  our  objective  is  simply  to  estimate  the  range  of  error  in  determining 
a  .(t),  Eq.  (10)  is  believed  to  be  adequate  because  Rocket  A31.  603  spins  and  tumbles 
only,  almost.  As  pointed  out  in  the  previous  paragraph,  Aa.(t)  is  expected  to  have 
its  mean  near  zero,  and  fora  ^(t),  Aa  ^(t)  varies  between  ±  !8°,  with  a  standard 
deviation  probably  lying  within  ±  10®  for  its  sinusoidal  behavior. 


7.  TEE  PROBLEM  OF  ATTITUDE  DETERMINATION  OF 
ROCKET  AXIS 

In  this  section,  we  determine  the  attitude  of  the  rocket  axis  with  respect  to  the 
local  vertical  coordinate  system.  For  Rocket  A31.  603,  magnetometer  No.  1  looks 
forward  along  the  spin  axis,  which  is  approximately  the  same  as  the  geometrical 
axis  of  the  rocket.  If  the  spin  axis  deviates  considerably  from  the  geometrical  axis, 
or  if  magnetometer  No.  1  dons  not  look  along  either  axis,  there  would  be  a  con¬ 
siderable  source  of  error  in  the  analysis.  Nevertheless,  even  as  a  zeroth  order 
approximation,  let  us  try  to  derive  as  much  information  about  the  rocket  attitude 
as  possible,  in  order  to  gain  an  understanding  of  the  behavior  of  the  rocket. 
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Let  us  summarize  what  we  know  at  this  stage.  By  using  two  magnetometer 
data  sets,  the  aspect  angle  ot . (t)  between  a  magnetometer  vector  and  the  Earth's 
magnetic  field  vector  D  at  the  poe  'tlon  of  the  rocket  is  known  approximately,  at 
Irregularly  spaced  instants,  about  1  sec  apart.  The  magnetic  field  vector  B  is  not 
sensitive  to  time  and  space  variation  for  the  duration  and  trajectory  span  of  the 
rocket  considered,  and  is,  therefore,  known  with  good  accuracy.  The  problem  now 
is  whether  one  can  determine  the  attitude  (azimuth  and  elevation  angles)  of  the  rocket 
axis  R  with  respect  to  a  local  vertical  coordinate  system  (x  =  east,  y  *  north, 
z  *  local  vertical,  that  is,  zenith). 


8.  ZEROTH  ORDER  APPROXIMATION  i 

Observation  of  a(t)  variation  against  t  shows  that  the  aspect  angle  of  the  rocket 
axis  varies  between  nearly  -180*  to  nearly  +180*,  implying  that  the  rocket  tumbles 
almost  in  a  vertical  plane.  The  deviation  from  the  plane  of  motion  is  approximately 
varying  at  less  than  15®.  That  is,  the  rocket  precesses  at  a  large  cone  angle, 
greater  than  75®,  so  that  it  behaves  almost  like  tumbling.  From  energy  considera¬ 
tions,  it  is  more  favorable  for  a  rocket  to  tumble  in  its  trajectory  plane  than  to 
wobble  significantly  out  of  its  trajectory  plane,  A  somewhat  analogous  situation  is 
known  in  the  theory  of  satellite  boom  dynamics,  in  which  modes  of  boom  oseilla- 

E  O 

tions  in  the  equatorial  plane  ’  of  the  dominant  angular  momentum  of  the  system 
are  energetically  more  favorable  than  those  wobbling  out  of  the  equatorial  plane. 

As  a  zeroth  order  approximation,  let  us  assume,  for  the  time  being,  that  the 
rocket  tumbles  in  its  vertical  trajectory  plane,  i’his  assumption,  together  with  the 
knowledge  of  a  j(t).  enables  us  to  constrain  the  possible  direction  of  the  rocket  axis 
from  those  lying  anywhere  in  a  cone  to  only  two  directions,  the  intersections  of  the 
a  -cone  with  the  vertical  plane.  Now,  we  have  obtained  two  sets  of  solutions,  at 
all  times  t. 

To  decide  which  solution  is  the  correct  one  at  time  t,  we  require  that:  (1)  the 


solution  can  be  traced  back  continuously  to  the  initial  time  tQ  so  that  the  two  sets  of 

solutions  can  be  distinctly  identified,  and  (2)  the  set  that  does  satisfy  the  initial 

condition  is  the  correct  one.  The  initial  condition  is  at  time  t  ,  when  the  rocket 

o 

starts  to  tumble,  due  to  the  separation  of  its  booster.  At  that  moment,  the  rocket 
is  still  pointing  along  its  velocity  vector.  For  Rocket  A31.  603,  we  have  found  that 
one  of  the  two  sets  of  solutions  satisfies  the  initial  condition  exactly.  This 


4 

* 

i 

* 


3  1 


5.  Lai,  S.T.,  ai ..£■  Bkavnani,  K.  H.  (1975)  Dynamics  of  Satellite  Wire-Boom  Systems, 
AFCRL-TR-75-U230,  AD  A014859.  - 

8,  Lai,  S.T.,  Mahon,  K. ,  and  Smiddy,  M.  (1979)  Dynamics  of  wire  boom  oscilla¬ 
tions  on  a  spinning  satellite,  J,  Appl.  Math.  Phys.  30:1-29. 
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encouraging  result  leads  one  to  have  confidence  on  the  correctness  of  the  selection 
of  the  solution  sets,  and  to  carry  on  further  for  refinement  of  results. 

0.  ROCKET  ELEVATION  AND  AZIMUTH  ANGLES 

We  define  the  azimuth  0  and  elevation  0  angles  in  a  local  vertical  coordinate 
system  as  in  Figure  1. 


Figure  1.  Local  Vertical 
Coordinate  System 


OR  is  the  vector  of  the  rocket  axis.  In  vector  notation,  we  write  the  unit  vector: 

* 

R  =  (cos  0  sin  p.  cos  0  cos  0.  sin  0)  .  <  1 1' 

A 

Similarly,  the  magnetic  field  vector  B  can  be  written  as 

A 

B  *  (cos  l  sin  O,  cos  l  cos  D,  sin  1)  (12) 

where  1  is  the  inclination  angle  and  D  the  dip  angle.  The  cap  *  denotes  unit  vector. 

At  White  Sands  where  the  rocket  considered  was  launched,  l  =  60*  and  D  =  11°. 

Also,  at  launch  time,  it  is  known  that  0  =  33.  5*  eastward  from  north.  It  is  believed 
that  this  angle  should  remain  unchanged  until  the  booster  separates  at  tQ  =  64  seconds. 

Since  the  rocket  tumbles  in  an  approximately  vertical  plane,  the  variation 
Max|0  -  <  0  >  I  should  be  small  compared  to  the  variation  of  &  in  a  tumble  period. 

For  the  rocket  considered,  the  ratio  of  such  variations  is  less  than  10%.  Thus, 
the  angle  0  behaves  almost  like  a  constant,  while  the  angle  0  varies. 

* 

From  Eqa.  (11)  and  (12),  the  angle  a  between  the  magnetic  field  vector  B  and 

A 

the  rocket  axis  vector  R  is  given  by 


. . . 


-  ■- !  - is!!  /i BSP?.'.  »?3pr.TJf7!!nV ' 


i-* 


r 


cos  <*  ■  B  *  3.  “  cos  0  cos  I  cos  (^  -  i»  +  sin  0  sin  I . 

a  % 

The  closest  encounter  between  B  and  R  occurs  when 
a  cos  o 


TT 


*  0 


(13) 


(14) 


which  gives  the  condition: 

•  .« 

.  tan  I 

tan  J  '  cSsiy  -~PT  * 


Substituting  Eq.  (15)  into  Eq.  (13),  we  have 


cos  <*min  *  |  cos2!  cos2  (^  -  D)  +  sin2l 


(15) 


r 


and 


^=Oi  cos 


-([ 


cos2amin-sin21 


j'7— } . 


(16) 


In  Eq.  (16).  there  are  two  sets  of  solutions.  We  can  choose  the  set  that  varies 
around  the  initial  value  at  t  *  t  where  $  *  33. 5*.  Thus,  this  technique  gives  the 
values  of  0  at  the  moments  of  closest  encounter  between  B  and  R,  that  is.  when 
a  1  is  minimum. 

For  the  rocket  considered,  there  are  only  six  moments  of  B  and  R  closest 
encounter.  Together  with  the  initial  value  of  tQ,  we  have  seven  data  points  to  be 
fitted  by  a  sine  curve  with  tumble  frequency.  The  fitted  function  gives  $(t)  at  all 
times.  Then,  using  Eq.  (13),  the  angle  0(t)  can  be  calculated  at  all  points.  This 
determines  the  elevation  and  azimuth  angles  of  the  rocket  at  all  points  (Figure  2). 
The  results  are  listed  in  Appendix  A. 

In  practice  there  i3  noise  in  the  data,  which  should  be  smoothed  by  curve-fitting 
or  other  filtering  techniques  before  the  method  of  closest  encounter  can  be  applied. 


10.  MOON  DIRECTION 

To  calculate  the  elevation  (e )  and  the  azimuth  (X)  of  the  moon  i 1  the  time  of 
rocket  flight,  we  use  the  data  in  Table  1. 


i,t»  dkik  ■ueatMUaMi 


Table  l.  Ephemeris  Data  of  Moon  and  Launch  Site 


Moon 

Greenwich 

White  Sands 

Right  Ascension 
Declination 

5h  54m  04* 

18. 47  5  • 

I7hllm  11* 

10h5m  11* 
32. 5» 

The  right  ascension  nl«<  ShS4mQ4*  -  •».  8*.  and  11*  *  Ml. »*•  *»  » 

coordinate  syntem,  E.  la  which  x  *  south.  y  •  east(both  on  Um  equatorial  piano), 
and  i  «  north  star.  th#  direction  cooinoa  of  aonith  i  rector  ia 

(  -  (coe  32. 8*.  0.  sin  32.8*) 

and  that  of  moon  M  rector  is 

j|  a  (cos  8  cos  ifa  cos  8  sin  tff,  sin  8 ) 

where 

if,  m  -l0h5mll*  +  5h54m04*  =  -82.8* 

8  -  18.475*  . 

In  a  local  vertical  system.  L.  in  which  x  =  south,  y  =  east,  and  s  =  senith,  e.  vector 
defined  in  E  system  becomes  matrix  transformed  by  a  rotation  about  -y-axts  of 
L -system.  That  is 

l  Ml,  ■  sin  32. 5*  0  -cos  32.5*  cos  8  cos  tf/ 

L  0  10  cos  8  sin  ^ 

cos  32. 5*  0  sin  32. 5*  sin  6 

-  I" -0.  034 

-0. 843 
-0. 536 

The  elevation  angle  (€ ).  measured  from  the  local  horison,  and  the  asimuth  angle 
(  X).  measured  clockwise  from  the  east,  are  then  obtained  by  equating 

M  ■  cos  t  sin  X 
cos  €  cos  X 
sin  e 

which  yields 

{«  =  32*4M  .  U7) 

(  X  »  182.3*  f 
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The  result*  in  Kq.  (1?)  hm  been  confirmed,  within  a  half  degree,  by  using  APGL's 
standard  computer  program  SOLUM.  They  do  net.  however,  agree  with  the  numbers 
given  by  Chamberlain7  and  by  Sluder  and  Kofsky. 8 

Using  the  rocket  attitude  data  and  the  moon  direction  [Eq.  (17)),  the  angle  0( t> 
between  the  rocket  axis  R  and  the  moon  g  at  time  t  can  be  calculated  (Figure  9). 
According  to  design,  the  moon  can  not  be  seen  unless  £(t)  lies  between  30*  and  44* 
(lor  backward  view)  or  between  176*  and  100*  (for  forward  view)  (see  Figure  4), 


The  range  of  moon  viewing  p  values  should  be  considerably  wider  than  those  spe¬ 
cified  by  design.  The  uncertainties  are  due  tot  (I)  a  angle  as  computed  by  using 
two  magnetometers.  (2)  spin  axis  deviation  from  rocket  axis.  (3)  noise  in  data, 

(4)  approximation  in  the  theory  of  B  and  R  close  encounter,  and  (5)  too  few  points 
available  in  the  close  encounter  curve  fitting.  The  uncertainty  due  to  source  (I) 
has  been  estimated  in  Section  8.  Let  us  put  up  a  rough  figure  for  the  combined 
uncertainty:  t  20*.  as  an  educated  guess  only,  unless  someone  can  give  a  reason 
to  indicate  otherwise. 

7 

The  moon  view  times  are  taken  from  Chamberlain  and  are  given  in  cc‘  nn  1 
of  Table  2.  The  magnetic  aspect  angles  oHt)  are  given  in  column  2.  The  ..nputed 
elevation  and  asimuth  angles  of  the  rocket  in  a  (east,  north,  senlth)  coordinate 
system  are  given  in  columns  3  and  4.  The  angles  f?  (t)  between  the  moonshine  vector 
and  rocket  vector  are  given  in  column  5. 


Table  2.  Rocket  Attitude  Angles  During  Moon  View  Times 


Moon  View 
Time,  sec 

Magnetic  Aspect 
Angle  a  (t) 

Rocket 

Elevation 

Rocket 

Asimuth 

Angle  p  (t) 
Between  Rocket 
and  Moon 

125.  1 

85.6* 

-125.5* 

31.  1* 

159.5 

38.8° 

-  39.3* 

58.  1* 

266.2 

36.7* 

-  97.5* 

43.  1* 

406.3 

21.2* 

-  75.3* 

52.9* 

46.3* 

Since  the  moonshine  is  almost  in  west  to  east  direction,  dipping  downwards,  it  is 
reasonable  that  the  rocket  should  lie  in  a  south-west  to  north-east  direction,  dipping 
downwards,  in  order  for  the  camera  on  the  rocket  to  view  the  moon.  Qualitatively, 
the  rocket  elevation  and  asimuth  angles  as  shown  in  Table  2  satisfy  this  criterion. 

As  for  the  angle  p  between  the  rocket  and  the  moon,  the  last  three  numbers  fall  is 
the  range  (30°,  44°)  ±  20*.  as  they  should.  The  first  j9  angle  lies  just  o*;tside  this 
range.  In  no  time  is  the  moon  seen  in  a  forward  view. 
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Appendix  A 

Attitude  Data  of  ItockM  No.  A31JQ3 

Column  1  showa  universal  time  starting  from  the  moment  of  booster  separation. 
Column  2  lists  magnetic  pitch  angles  [ref.  2).  Columns  3  and  4  are  calculated  ele¬ 
vation  and  asimuth  angles  in  the  local  vertical  coordinate  system.  Column  5  gives 
the  calculated  moon  pitch  angles. 
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4  4.05 

61.84 

*57 

*39.40 

*  *  .  4  P 

-  r  4 , 1  * 

45.4  5 

6  0.4  -1 

35  P 

733.80 

72.  24 

- r2 . 61 

45.  »  4 

59.  54 

360 

*33. 4C 

71.07 

-ci .71 

46.42 

50.45 

36  0 

400,00 

*0.0* 

- c- 0 .  0* 

46,3° 

6  7.5’ 

261 

4  *  1 . 8  C 

2  0  .  C  0 

-67.46 

4  >  .7  4 

5  5.6* 

362 

4"1.60 

26.  10 

-»4  ,  P5 

48.47 

57.74 

36  3 

402.40 

25,  0  3 

-83.28 

49.1  8 

52.5* 

36  4 

402. 2C 

24.  94 

-(3.05 

49.9  7 

52.31 

365 
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4 CM..  80 
406.8  C 

407.40 
408.00 

409.60 

400.60 

410.60 

411.20 

412.40 

413.40 
41 4 . 2  C 

414.60 

414.80 

416.60 
416.00 

416. 80 
41*. 20 

418.40 

418.80 
410.8 C 

420.80 
421. 4C 
422.00 
4*2.60 
423. 8C 

426.60 

427.40 
428.00 

423.60 

429.40 

426.80 
4*0. 60 

431.60 

438.40 

433.60 

434.80 

435.80 

436.20 
436.6  0 
437.00 

437.40 
438.00 

439.80 

439.20 
4  39.8  0 

440.40 

440.80 

441.80 

442.60 

443.80 
445.00 

446.80 

448.40 

448.80 

449.40 


23.08 

21.26 

21.05 

21.2? 

22.86 

23.27 

23.50 

23.93 

25.11 

26.  07 

26.75 
26.96 

27.  07 

28.  42 
29.07 
30.04 

30.72 
33.  CP 
*3.55 
35.14 

36.71 
37.  84 
3«?.  C5 
41.03 
43.  58 
51.59 

53.72 
54.  42 
57.10 

59.56 

60.90 
63.  55 

66 .77 
69.  3u 
73.02 
76,69 
79.  53 
80. 6* 

81.78 

82.91 
83.  8° 
85.38 

87.40 

88.37 
89.77 
91.18 
92.1? 
94.  27 
95.84 

97,54 

ICO, 11 

104.72 
109.15 
110. 28 
112.41 


-79.80 
*75.47 
-74.64 
-59.08 
-65.63 
-55.45 
-124.40 
-56.06 
-53.40 
-52.11 
-51.26 
-51.83 
-*0.92 
-49, 00 
-48.11 
-46. 85 
-45.94 
-42.46 
-42.26 
-40.27 
-36 , 31 
-*6.°0 
-35.  41 
-32.98 
-29. e9 
-20. 42 
-17.9* 
-16.94 
-13. 98 
-11.14 
-0.68 
-6.  55 
-2.85 
.04 
4,  31 
8.51 

11.78 
13.07 
14.  37 
15.67 
16.  81 

18.54 
20.87 
22.  00 
23. 63 

25.27 

26.37 
*8.  86 
70.7? 
72,81 
76  •  7° 

41.12 
4*.  99 
4*. *6 
49.59 


A7I9UTH 

066, 

51.20 
52.73 

53.16 

53.57 
53.97 

54.61 

55.20 

55.54 

56.16 

56.5  2 
56.96 

57.12 
57.19 
5  7.46 

67.61 
57.84 
57.99 
56.22 

68.29 
56.44 

56.54 
5A.57 
5  8.5  8 

58.57 

56.5  1 
58,07 
57.87 
57.70 

57.51 
57.23 
57.06 

56.75 

56.30 
55.9  1 

56.26 

54.55 

53.91 
53.64 

53.37 
53.09 

52.80 
52.36 

51.75 
51.44 
50.96 

50.46 
50. 1  3 
49.»7 

48,56 

47.46 
46,32 

44.5  6 
42.93 

42.52 
41.90 


R.MOCN 

rec. 


fom 

NO. 


49.82 

366 

46.49 

367 

45.8? 

*68 

35.94 

369 

34.15 

370 

33.65 

371 

87.02 

372 

33.00 

37  3 

31.90 

374 

31.07 

*75 

3C.51 

376 

3C.  33 

377 

3  0.24 

378 

*9.  30 

379 

2B.90 

380 

2?  •  *? 

381 

27.97 

38  2 

27.  0  2 

38* 

26.83 

354 

26.41 

385 

26.  16 

*86 

26.IJQ 

*87 

26.11 

38  8 

26.34 

386 

26.96 

29  C 

30.88 

391 

32.34 

39? 

3  3.59 

39* 

34.91 

39  4 

?  6.  91 

3S5 

*8.04 

39  6 

40,36 

39  7 

4  3.32 

39e 

45.73 

?96 

49.41 

400 

53.1? 

401 

56.1? 

40  2 

57. *1 

403 

58.  51 

40  4 

59,72 

406 

60.81 

406 

62.46 

40  7 

6  4.69 

4oe 

65.78 

409 

67.38 

410 

68.  97 

411 

7C.04 

412 

72.5* 

413 

74.40 

414 

76.63 

415 

79.60 

416 

eu.6* 

417 

89.27 

418 

90.44 

419 

92.52 

420 
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